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\Q ■ Abstract. We consider the double scaling limit in the random matrix ensemble with an 

external source 

O . 1 e -nTr(lM*-AM) dM 

defined on n x n Hermitian matrices, where A is a diagonal matrix with two eigenvalues ±a 
D ' of equal multiplicities. The value a = 1 is critical since the eigenvalues of M accumulate as 

n — > oo on two intervals for a > 1 and on one interval for < a < 1. These two cases were 
treated in Parts I and II, where we showed that the local eigenvalue correlations have the 
universal limiting behavior known from unitary random matrix ensembles. For the critical 
case a = 1 new limiting behavior occurs which is described in terms of Pearcey integrals, as 
shown by Brezin and Hikami, and Tracy and Widom. We establish this result by applying 
t^J- ■ the Deift/Zhou steepest descent method to a 3 x 3-matrix valued Riemann-Hilbert problem 

, which involves the construction of a local parametrix out of Pearcey integrals. We resolve 

the main technical issue of matching the local Pearcey parametrix with a global outside 

, parametrix by modifying an underlying Riemann surface. 

\o : 
o ■ 

p ,1 l. Introduction and statement of results 

1.1. The random matrix model. This is the third and final part of a sequence of papers 
on the Gaussian random matrix ensemble with external source 

>• z n K J 

^ ■ defined onnxn Hermitian matrices, where A is a diagonal matrix with two eigenvalues ±a 
(with a > 0) of equal multiplicities (so that, n is even). This matrix ensemble was introduced 
by Brezin and Hikami [5J CD] as a simple model for a phase transition that is expected to 
exhibit universality properties. The phase transition can be seen from the behavior of the 
eigenvalues of M in the large n limit, since for a > 1, the eigenvalues accumulate on two 
intervals, while for < a < 1, the eigenvalues accumulate on one interval. The limiting mean 
density of eigenvalues follows from earlier work of Pastur [22] • It is based on an analysis of 
the equation (Pastur equation) 

e-ze + (l-a 2 )£ + a 2 z = 0, (1.2) 

which yields an algebraic function £(z) defined on a three-sheeted Riemann surface. The 
restrictions of £(z) to the three sheets are denoted by £j(z), j = 1,2,3. There are four real 
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Figure 1. The structure of the Riemann surface for the equation (jl.2j) for 
the values a > 1 (left), a = 1 (middle) and a < 1 (right). In all cases the 
eigenvalues of M accumulate on the interval (s) of the first sheet with a density 
given by (jl.3|) . 



branch points if a > 1 which determine two real intervals. The two intervals come together 
for a — 1, and for < a < 1, there are two real branch points, and two purely imaginary 
branch points. Figure [T] depicts the structure of the Riemann surface £(z) for a > 1, a = 1, 
and a < 1. 

In all cases we have that the limiting mean eigenvalue density p(x) = p(x; a) of the matrix 
M from (jl.lj) is given by 

p(x; a) = — Im£i + (x), x £ R, (1.3) 

71 

where £,i+(x) denotes the limiting value of £1(2) as z — > x with Imz > 0. For a = 1 the 
limiting mean eigenvalue density vanishes at x = and p(x; a) ~ Ix] 1 ^ 3 as x — > 0. 

We note that this behavior at the closing (or opening) of a gap is markedly different from 
the behavior that occurs in the usual unitary random matrix ensembles Z~ l e - nTrV ( M ) dM 
where a closing of the gap in the spectrum typically leads to a limiting mean eigenvalue 
density p that satisfies p(x) ~ (x — x*) 2 as x — > x* if the gap closes at x = x*. In that case 
the local eigenvalue correlations can be described in terms of ^-functions associated with the 
Painleve II equation, see The phase transition for the model (jl.lj) is different, and it 

cannot be realized in a unitary random matrix ensemble. 

1.2. Non-intersecting Brownian motion. The nature of the phase transition at a = 1 
may also be seen from an equivalent model of non- intersecting Brownian paths, see Figure El 
Consider n independent one-dimensional Brownian motions that are conditioned to start at 
time t — at the origin, end at time t — 1 at ±1, where half of the paths ends at +1 and the 
other half at —1, and that are conditioned not to intersect at intermediate times t G (0, 1). 
As explained in [2], at any intermediate time t, the positions of the n Brownian motions, 
have the same distribution as the eigenvalues of a Gaussian random matrix ensemble with 
external source (up to trivial scaling). Now, as n —>■ oo and under appropriate scaling of 
the variance of the Brownian motions, the paths fill out a region in the t — x-plane. Then 
for small time the paths are in one group, which at a certain critical time t cr splits into two 
groups, where one group ends at x — +1 and the other group at x = —1. The situations 
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t cr , and t > t c 



correspond to a < 1, a = 1, and a > 1, respectively, in the 
Gaussian random matrix model with external source. 

The boundary curve has a cusp singularity at the critical time as shown in Figure El 




Figure 2. Non- intersecting Brownian paths that start at one point and end at 
two points. At any intermediate time the positions of the paths are distributed 
as the eigenvalues of a Gaussian random matrix ensemble with external source. 
As their number increases the paths fill out a region whose boundary has a 
cusp. 



1.3. Correlation kernel. Brezin and Hikami jUJ EI! showed, see also [2E] 5 that the eigen- 
values of the random matrix ensemble (jl.lj) are distributed according to a determinantal 
point process. There is a kernel K n (x, y \ a) so that the eigenvalues xi, . . . , x n have the joint 
probability density 

p n (xi, ...,x n ) = — det(K n (xj , x k ; a)) jtk =i,...,n 

and so that for each m < n, the m-point correlation function 

n\ f f 

Rmi^Xi • • • ? •^m) 7 \T I I Pni-^lj ■ ■ ■ j •^n)dx m j r \ • • • dx n 

[n-m)\ J J 

n—m times 

takes determinantal form as well: 

Rm{%li • • • i %m) det (-K"n [Xj , Xfc, (3.) ) j' J fc=l,...,m- 

In jH] we pointed out that the kernel can be built out of multiple Hermite polynomials 
|2l El I2S] in much the same way that the correlation kernel for unitary random matrix 
ensembles (without external source) is related to orthogonal polynomials. The Christoffel- 
Darboux formula for multiple orthogonal polynomials [HI E] allows one to express the kernel 
in terms of the Riemann-Hilbert problem for multiple Hermite polynomials (see j2Z] and 
below). Applying the Deift/Zhou steepest descent analysis ITH] to the Riemann-Hilbert 
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Figure 3. The contour S that appears in the definition of q(y). 



problem in the non-critical case, we were able to show that the kernel has the usual scaling 
limits from random matrix theory. That is, we obtain the sine kernel 



K^(x,y) = Sini ; {X -y ) (1.4) 
n(x - y) 

in the bulk, and the Airy kernel 

K ed 9 e^ s = Ai(x)Ai'(?/) - Ai' (x)M(y) 

x-y 

at the edge of the spectrum, as scaling limits of K n (x, y; a) if a > 1 j7j or < a < 1 0. 

1.4. Double scaling limit. In this paper we consider the double scaling limit at the critical 
parameter a = 1 of the Gaussian random matrix ensemble with external source, or equiva- 
lently, of the non-intersecting Brownian motion model at the critical time t = t cr . As is usual 
in a critical case, there is a family of limiting kernels that arise when a changes with n and 
a ^ 1 as n ^ oo in a critical way. These kernels are constructed out of Pearcey integrals 
and therefore they are called Pearcey kernels. The Pearcey kernels were first described by 
Brezin and Hikami [§J EI] . A detailed proof of the following result was recently given by 
Tracy and Widom [25] . 

Theorem 1.1. We have for every fixed b G R, 

lim (-^jj, -L; 1 + ttV) = K cusp (x,y;b) (1.6) 



n ^oo n 3/4 " ^ n 3/4> n 3/4 ' 2 ^ 

where K cusp is the Pearcey kernel 



r^cuspf t \ P(x)q (y) ~ p (x)q (y) +p jx)q{y) - bp(x)q(y) 

K 1J [x,y,b) = (1.7) 

x-y 



with 



p ( x ) = — f e--^-^ s2+isx ds and q(y) = — f e^+^+^dt. (1.8) 

The contour E consists of the four rays argy = ±7r/4, ±37r/4, with the orientation shown in 
Figure El 
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The functions (|1.8jl are called Pearcey integrals j2Bj- They are solutions of the third order 
differential equations p"'(x) = xp(x) + bp'(x) and q"'(y) = —yq(y) + bq'(y), respectively. 
Away from the critical point x = 0, the usual scaling limits (jl.4|) and ()1.5|) from random 
matrix theory continue to hold in the case a = 1 (also in the double scaling regime). This 
can be proved for example as in 01213 , and we will not consider this any further here. 

Theorem 11.11 implies that local eigenvalue statistics of eigenvalues near are expressed in 
terms of the Pearcey kernel. For example we have the following corollary of Theorem ll.il 

Corollary 1.2. The probability that a matrix of the ensemble (jl.ip with a = 1 + bn~ 1 / 2 /2 
has no eigenvalues in the interval [cn~ 3 / 4 , dn~ 3 ^ 4 } converges, as n — ^ oo ; to the Fredholm 
determinant of the integral operator with kernel K cusp (x, y; b) acting on L 2 (c, d). 

Similar expressions hold for the probability to have one, two, three, . . . , eigenvalues in an 
0(n~ 3 / 4 ) neighborhood of x = 0. 

Tracy and Widom j2H] and Adler and van Moerbeke gave differential equations for 
the gap probabilities associated with the Pearcey kernel and with the more general Pearcey 
process which arises from considering the non-intersecting Brownian motion model at several 
times near the critical time. See also [20] where the Pearcey process appears in a combina- 
torial model on random partitions. 

I. 5. Steepest descent method for RH problems. Brezin and Hikami and also Tracy 
and Widom used a double integral representation for the kernel in order to establish Theorem 

II. 11 In this paper we use the Deift/Zhou steepest descent method for the Riemann-Hilbert 
problem for multiple Hermite polynomials. This method is less direct than the steepest 
descent method for integrals. However, an approach based on the Riemann-Hilbert problem 
may be applicable to more general situations, where an integral representation is not avail- 
able. This is the case, for example, for the general (non-Gaussian) unitary random matrix 
ensemble with external source 

_L e -nTr(V(M)-AM) dM (Lg) 

Z n 

with a general potential V. The Riemann-Hilbert problem is formulated in Section 2. 

The asymptotic analysis of the Riemann-Hilbert problem presents a new feature that we 
feel is of importance in its own right. We will not use the Pastur equation (jl.2j) which 
defines the ^-functions and the Riemann surface that corresponds to it, but instead we use a 
modified equation to define the ^-functions. We discuss this in Section 3. The modification 
may be thought of in potential theoretic terms and we briefly discuss this in Section 3 as 
well. 

The anti-derivatives of the modified ^-functions are introduced in Section 4 and they play 
an important role in the steepest descent analysis of the Riemann-Hilbert problem in the 
rest of the paper. The main issue is the construction in Section 8 of the local parametrix 
around with the aid of Pearcey integrals. The modification of the ^-functions is used here 
to be able to match the local Pearcey parametrix with the outside parametrix. Even so it 
turns out that we cannot achieve the matching condition on a fixed circle around the origin, 
but only on circles with radii n~ 1//4 that decrease as n increases. However, the circles are big 
enough to capture the behavior (jl.6|) which takes place at a distance to the origin of order 
n -3/4 rj^ e p rec j se estimates that lead to the proof of Theorem 11.11 are given in the final 
Sections 9 and 10. 
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2. RlEMANN-HlLBERT PROBLEM 

As shown in our paper [B], the correlation kernel is expressed in terms of the solution to 
the following 3x3 matrix valued Riemann-Hilbert (RH) problem. 
Find Y : C \ E -> C 3x3 such that 

• Y is analytic on C \ 1, 

• for x G M, we have 





g—n(^x 2 —ax) 


e -n{\x 2 





1 











1 



Y+(x) = Y-(x)\0 1 (2.1) 

\0 I J 

where Y + (x) (Y_(x)) denotes the limit of Y(z) as z — > x from the upper (lower) 

half-plane, 

as z —>■ oo, we have 

K(z) = ( J + OI-) ) I I . (2.2) 




where / denotes the 3x3 identity matrix. 

The RH problem has a unique solution, given explicitly in terms of the multiple Hermite 
polynomials. The correlation kernel of the Gaussian random matrix model with external 
source is equal to 

Kn{x > m a) = 2m(x - y) (° ^ e ~ na ^ Y + M Y +W \ ° (2 - 3) 

In what follows we are going to apply the Deift/Zhou steepest descent method for RH 
problems to the above RH problem for Y. It consists of a sequence of explicit transformations 
Y i— > T i— > S i— > R which leads to a RH problem for R in which all jumps are close to the 
identity matrix and which is normalized at infinity. Then R is close to the identity matrix, 
and analyzing the effect of the transformations on the kernel ()2.3|) we will be able to prove 
Theorem 11.11 



3. Modification of the ^-functions 

3.1. Modified Pastur equation. The analysis in [21 [7] for the cases a > 1 and < a < 1 
was based on the equation (jl.2j) and it would be natural to use (jl.2j) also in the case a — 1. 
Indeed, that is what we tried to do, and we found that it works for a = 1, but in the double 
scaling regime a = 1 + with b ^ 0, it led to problems that we were unable to resolve in a 
satisfactory way. A crucial feature of our present approach is a modification of the equation 
(jl.2j) when a is close to 1, but different from 1. At x = we wish to have a double branch 
point for all values of a so that the structure of the Riemann surface is as in the middle 
figure of Figure d for all a. 

For c > 0, we consider the Riemann surface for the equation 
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where w is a new auxiliary variable. The Riemann surface has branch points at z* = ^rpc, 
—z* and a double branch point at 0. There are three inverse functions Wk, k = 1,2, 3, that 
behave as z — ► oo as 

c 2 / 1 

= Z hO — 

"•W= e + S + °(?) (3 ' 2) 

and which are defined and analytic on C \ [— z*, z*], C \ [0, z*] and C \ [—2*, 0], respectively. 
Then we define the modified ^-functions 

€k = Wk + —, for k = 1,2, 3, (3.3) 
which we also consider on their respective Riemann sheets. In what follows we take 

and p = c 2 -l. (3.4) 



a + a/o 2 + 8 , > 



Note that a = 1 corresponds to c = 1 and p = 0. In that case the functions coincide with 
the solutions of the equation (jl.2|) that we used in our earlier works. From (|3.1j) . ()3.3|) . and 
(|3.4|) we obtain the modified Pastur equation 

- zi 2 + (1 - a 2 )£ + a 2 z + ^~ 2 ^ = 0, (3-5) 

where c is given by ()3.4j) . 

Lemma 3.1. Lei a > and tafce c and p as in (|3.4|) . Tnen at infinity we have 

^ + o(i), (3.6) 

Proof. This follows from direct calculations using ()3.2j) . (j3.3|) and the fact that 2c — - = a. 
Alternatively, one could also use ()3.5|) . □ 

The new ^-functions have the same asymptotic behavior ()3.6|) as z — ► 00 (up to order 
1/z 2 ) as the solutions of ([1.2)1 . This is important for the first two transformations of the 
Riemann-Hilbert problem. The situation at z = is different. The fact that we can control 
the behavior at z = as well is the reason for the introduction of the modified ^-functions. 
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3.2. Behavior at z = 0. We start with the behavior of the functions Wk- 

Lemma 3.2. There exist analytic functions f\ and g\ defined in a neighborhood U\ of z = 
so that for z G U\ and k = 1, 2, 3, 

r _jik z i/z^ ^ _ u k z 5/3 gi ^ + £ f O rlmz>0, 

w k (z)=\ 3 (3.7) 

J -uj k z 1,3 f l (z) - uo 2k z 5/3 9l (z) + - for Im z < 0. 

In addition, we have /i(0) = c 2 / 3 , and fi(z) and g\{z) are real for real z G U\. 

Proof. Putting z = x 3 and w = xy in (|3.1|) we obtain 

y 3 = x 2 y 2 - c 2 (3.8) 

which has a solution y = y(x) that is analytic in a neighborhood U\ of and satisfies 
y(0) = — c 2 / 3 and y'(0) = 0- Then we can write y(x) = —fi(x 3 ) — x A gi(x 3 ) + x 2 hi{x 3 ) 
with fi, gi and hi analytic in U\ and /i(0) = c 2//3 . Putting this back in ()3.8j) we find after 
straightforward calculations that (with z = x 3 ) 

f 1 {z)g l {z)= 1 -, f 1 ( z ) 3 -c 2 + ^z 2 + g 1 (z) 3 z' = 0, (3.9) 

and hi(z) = |. Going back to z and w variables, we see that there is a solution w = w(z) to 
flUD with 

w(z) = -z^h{z) - z^ 3 9l (z) + |, for z E C \ (-oo, 0], 

where we take the principal branches of the fractional powers. This solution is real for z real 
and positive, and so it coincides with the solution w^{z). This proves ()3.7|) for k = 3. The 
expressions f!3.7|) for k — 1,2 follow by analytic continuation. 

Since y(x) is real for real x, we also find that fi{z) and gi(z) are real if z is real. □ 

From ()3.9|) it is easy to give explicit expressions for f\ and g\. However we will not use 
this in the future. 

From Lemma f3. 21 and ()3.3|) we get the following behavior for the functions near z = 0. 

Lemma 3.3. There exist analytic functions f2 and g<i defined in a neighborhood U2 of z = 
so that for z 6 U2 and k = 1,2, 3, 

-tu 2k z 1/3 f 2 (z) - cu k z- 1/3 g 2 (z) + - for Im z > 0, 

= { I (3-10) 

-uo k z l/3 f 2 {z) - uo 2k z- l/3 g 2 {z) + - /or Im 2 < 0. 

In addition, we have 

h(0) = c 2/3 + V 4 / 3 (c 2 - 1), 3,(0) = c- 2 l\c 2 - 1), (3.11) 

and f 2 (z) and g 2 {z) are real for real z G U 2 . 

Proof. This follows from ()3.3|) and the previous lemma. Indeed from ()3.7|) where f\ and g\ 
satisfy the relations (|3.9|) . we can deduce 

. z W±-(f x {z) + 3z 2 0l (z) 2 ) - z-wl-Qf^z) 2 + z 2 0l (z)). 



w 3 (^) 3c 2 3c 2 
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Then ()071jl follows from flH2|) and flEJJ if we take 

Mz) = h{z) + ^f^(fi(z) + 3z 2 9l (z) 2 ) (3.12) 

and 

<&(*) = z 2 9l (z) + ^(S/^) 2 + ^)). (3.13) 

Because of Lemma (3.21 this also implies ()3.11j) and the fact that f 2 {z) and #2(2) are real for 
real z G U 2 . □ 

3.3. Potential theoretic interpretation. As an aside we want to mention that the mod- 
ified ^-functions may be thought of in terms of a modified equilibrium problem for logarith- 
mic potentials. For a > 1, it was noted in that the limiting mean eigenvalue density 
p(x) = p(x; a) may be characterized as follows. We minimize 

E(px,p 2 ) = / / log] — - — rdfii(x)dfi 1 (y) + [ [ log- — - — :dp 2 (x)dp 2 (y) 
JJ \x-y\ JJ \x-y\ 

+ J J log j—^—^dp 1 (x)dp 2 (y) + J (^x 2 - ax^j dp^x) (3.14) 

+ J (j^x 2 + ax ^j dp 2 (x) 

among all non- negative measures pi, p 2 on M with J dpi = J dp 2 = ~. There is a unique 
minimizer [21|, and for a > 1, we have that supp(/ii) C [0, 00), supp(/i 2 ) C (— 00, 0], and p 
is the density of p\ + p 2 . For a < 1, the minimizing measures for (|3.14|) do not have disjoint 
supports, and in fact these minimizers are not related to our random matrix ensemble (jl.lj) 
at all. 

The modification we are alluding to is to minimize ()3.14|) among signed measures p\ = 
Pi — Pi , p 2 = p 2 — p 2 where pf are non- negative measures, that satisfy J dpi = J dp 2 = | 
and in addition 

(1) supp(/ii) C [0,oo), supp(/i 2 ) C (—oo,0], and, 

(2) There is a 5 > 0, such that supp(/i^) C [0,5) and supp(yU 2 ) C (—5,0]. 

The condition (1) plays a role for a < 1, since it prevents the supports of pi and p 2 to 
overlap. For a > 1, the condition (2) plays a role, since it allows the measures to become 
negative near 0. Now let pi, p 2 be the minimizers for this modified equilibrium problem, and 
let p be the density of p\+ p 2 . Then it can be shown that the density of p\ + p 2 is equal to 
^Im£i + (a;) where £1 is the modified £i-function introduced in this section. 

We will not use this potential-theoretic connection in the analysis that follows in this 
paper, but we anticipate that it might be important for the general unitary random matrix 
ensemble with external source (jl.9|) . 

We finally note that a modified equilibrium problem was also used in [iT] IT2"] in order 
to analyse the double scaling limit in unitary random matrix ensembles (without external 
source), so one might speculate that such an approach might be characteristic for double 
scaling limits in random matrix ensembles. 
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4. The A-functions 

4.1. Definition and first properties. The main role is played by the A-functions which 
are anti-derivatives of the ^-functions. They are defined here as 

h{z) = [ £k(s)ds (4.1) 

where the path of integration starts at on the upper side of the cut and is fully contained 
(except for the initial point) in C \ (—oo,z*] for k = 1,2, and in C \ (— oo,0] for k = 3. 
Then Ai and A2 are defined and analytic on C \ (—00, z*], and A3 is defined and analytic on 
C\(-oo,0]. 

As follows from ()3.6|) and (|4.1|) . the A-functions behave at infinity as 

\ 1 (z) = ±z 2 -\ogz + £ 1 + 0(l/z), 

X 2 (z) = az + Uogz + £ 2 + 0(1/ z), (4.2) 
X 3 (z) = -az + ^ \ogz + £ 3 + 0(1/ z), 

for certain constants £k, k — 1, 2, 3, where logz is taken as the principal value, that is, with 
a cut along the negative real axis. 

From contour integration based on ()3.6|) where we use the residue of £2 at infinity, we find 
Ai-(0) = ni and A2-(0) = —ni. Then we get the following jump properties of the A-functions 
on the cuts (— 00, 0] and (—00,2*]: 

Ai+ = A 2 - + ni, A 2 + = Ai_ - ni, A 3+ = A 3 _ on [0, z*], 

A 1+ = A 3 _, A 2+ = A 2 _ + ni, A 3+ = Ai_ - ni, on [-z*, 0], (4.3) 

Ai + = Ai_ — 2ni, A 2+ = A 2 _ + ni, A 3+ = A 3 _ + ni, on (-00, -z*}. 

4.2. Behavior near z = 0. Near the origin the A-functions behave as follows. 

Lemma 4.1. There exist analytic functions / 3 and g 3 in a neighborhood U3 of z = so that 

-^u 2k z^f 3 (z) - \uj k z^g,( z ) + Z 1 forlmz> 0, 

Xk _( ) - -uj k z^f 3 (z) - -uj 2k z^g 3 {z) + - for \mz < 0, 
In addition, we have 

h(0) = /a(0) = c 2 / 3 + l -c-V\c 2 - 1), g 3 (0) = 3g 2 (0) = 3c^ 3 (c 2 - 1), (4.5) 
and f 3 (z) and g 3 {z) are real for real z EU 3 . 

Proof. The relations ()4.4j) follow by integrating ()3.1()j) . Note that Ai_(0) = ni, A2-(0) = —ni, 
and A 3 _ (0) = 0. The other statements of the lemma also follow directly from Lemma 1331 □ 
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Rep.,) = Re(l 2 ) 

Rell,) = Re(l 3 ) 

_ Re(X 2 ) = Re(l 3 ) 



Figure 4. Curves where Re A., = ReA^ for the value a = 2.0. 

4.3. Critical trajectories. Curves where Re A., = ReA/c for some j ^ k are shown in 
Figures IH 03 and|HJ for the cases a > 1, a = 1, and a < 1, respectively. These are critical 
curves that play a crucial role in the asymptotic analysis. The curves are critical trajectories 
of the quadratic differentials (£j(z) — ^k{z)) 2 dz 2 (and their analytic continuations beyond the 
branch cuts in case a < 1). 

The solid curves in Figures EHH1 are the critical trajectories of the quadratic differential 
(£1(2) — £,2(z)) 2 dz 2 . The quadratic differential has a simple zero at z = z*. Three trajectories 
are emanating from z = z* at equal angles, one of these being the real interval (0, z*). For 
a > 1, the quadratic differential has a double zero at z = x for some x G (0,z*). Four 
trajectories are emanating from the double zero at equal angles as can be seen in Figure 

The dashed curves are the critical trajectories of the quadratic differential (£i(z)— ^{z)) 2 dz 2 . 
Because of symmetry, these are the mirror images of the trajectories of the quadratic differ- 
ential (£1(2) — (,2(z)) 2 dz 2 with respect to the imaginary axis. For a > 1, the solid curve that 
passes vertically through x an d its dashed mirror image with respect to the imaginary axis 
meet in two points ±iy on the imaginary axis. Together they enclose a neighborhood of the 
origin. 

The dashed-dotted curves are the critical trajectories of (£2(2) — £,z{z)) 2 dz 2 . For a < 1, the 
quadratric differential has two double zeros at z = ±iyo for some yo > 0. Four trajectories 
are emanating from these double zeros at equal angles as shown in the Figure |H1 Besides 
the imaginary axis there are curves passing horizontally through ±iyo, and these curves 
meet each other at two points ±Xo on the real axis and they enclose a neighborhood of 
the origin. Beyond these two points the quadratic differentials have analytic continuations, 
but the formula changes since either £ 2 or £ 3 reaches its branch cut and changes into £1. 
Consequently, the dashed-dotted curves in Figure El continue beyond ±x as either solid or 
dashed curves. 

The relative orderings of the real parts ReAi, ReA 2 and Re A3 changes if we cross one 
of the critical trajectories, but it remains constant in the regions bounded by the critical 
trajectories. For each of the unbounded regions we can determine the ordering from the 
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Re(V,) = Re(l 2 ) 

He(XJ = Re(l 3 ) 

_ Re(y = Re(l 3 ) 



Figure 5. Curves where Re A., = ReA^ for the value a = 1.0. 

behavior at infinity (|4.2jl . For example, we have in the right-most region ReAi > ReA2 > 
ReA 3 , and if we cross the solid curve where ReAi = ReA 2 , the ordering becomes ReA 2 > 
Re Ai > Re A 3 , and so on. 

In the cases a < 1 and a > 1 the trajectories enclose a bounded neighborhood of the 
origin. There is no such neighborhood in case a = 1. The neighborhood is small if a is close 
to 1. In this neighborhood the relative ordering of the real parts is different. 

So we can easily verify the following. 

Lemma 4.2. Except for z in the exceptional bounded neighborhood of the origin, we have 
that 

ReA 2 (-2) > max(Re Ai(^), Re As(z)) 
in the region in the right-half plane, bounded by the solid and dashed-dotted curves, and 

ReA 3 (z) > max(Re Xi(z), Re A 2 (z)) 

in the region in the left- half plane bounded by the dashed and dashed-dotted curves. 

The exceptional neighborhood will not cause a problem to us, since it turns out to shrink 
fast enough if a = 1 + {b/2)n~ l l 2 and n — > oo. For n large enough, the exceptional neighbor- 
hood is well within the disk around the origin of radius n -1 / 4 where we are going to construct 
a special parametrix with Pearcey integrals. Then the different ordering of the real parts of 
the Xk will not play a role. 

5. First two transformations of the RH problem 

The first and second transformation of the RH problem are the same as in our earlier 
paper [7j, except that we use the A-functions that were introduced in the last section via the 
modified ^-functions. 
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Figure 6. Curves where Re A., = ReA^ for the value a = 0.5. 

5.1. First transformation Y i— > T. Using the functions A& and the constants £k, k = 1, 2, 3, 
we define 

T(» = diag (e- nh ,e- nl \e~ nh ) F(z)diag (e^M*)- 1 ^ 2 ) ^ e n{x 2 {z)-az) ^ e n(x 3 (z)+az)\ _ ^ 

Then by (12. lj) and (|5.1|) and the jump properties f)4.3|) we have T + (x) = T_(x)j T (x) for 
x G R, where 

( g n(Ai-A 2 ) + ^ e n(A 3 -Ai_) N 

e n(A!-A a )_ o J 16(0,4 (5.2) 

1 

( e n(Ai-A 3 )+ e n(A 2 +-Ai_) ^ 
1 J. ,■£(_,-. 0). (5.3) 

e n(Al " A3) - 

(]_ e n(A 2+ -Ai_) e n(A 3+ -Ai_)\ 
1 x G (— oo, —-2*) U (z* , oo). (5.4) 

1 J 

The function T(z) solves the following RH problem: 

• T is analytic on C \ R, 

• for x G R, we have 

T + (x) = T_(x)j T (x), (5.5) 
where jj* is given by (|5.2|) - ([5.4|h 



• as z — > cx), 

r(z) = I + of~Y (5.6) 

The asymptotic property ()5.6|) follows from ()2.2j) . (|5.1|) . and the behavior ()4.2|) of the A- 
functions at infinity. 
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Figure 7. Opening of lenses around the intervals [0,z*] and [— z*,0] for the 
value a = 2.0. The upper and lower lips of the lenses together with the real 
axis form the contour £5 which is shown in bold. Also shown are the critical 
trajectories as in Figure |U 



5.2. Second transformation T 1— > S. The second transformation of the RH problem con- 
sists of opening of lenses around the intervals [0, z*\ and [— z*,0]. The lenses are as shown 
in Figure H We define (see also Section 5 in [7j) 



1 

S = T ( - e n ( Al " A2 ) 1 - e n ( A 3- A 2) J in the upper right lens region, (5.7) 

1 

1 

S = T ( e n(Al ~ A2) 1 -e n ( A3 ~ A2 ) I in the lower right lens region, (5.8) 

1 

1 N 
S — T I 1 ] in the upper left lens region, (5.9) 

_ e n(Ai— A3) _ e n(A 2 -A 3 ) j 

1 N 

S — T I 1 ] in the lower left lens region, (5.10) 

3 n(Ai-A 3 ) _ e n(A 2 -A 3 ) ^ 



and S = T outside the lenses. 

It leads to a matrix valued function S which is defined and analytic in C \ £5, where 
consists of the real line and the upper and lower lips of the lenses. On E5 we have S + = S-js 
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where js is defined as follows (the orientation on Eg is taken from left to right): 

on (0,2*), (5.11) 



on(-z*,0), (5.12) 



on (-oo, -z*) U (**,oo), (5.13) 



on the upper lip of the right lens, (5.14) 



on the upper lip of the left lens, (5.15) 



on the lower lip of the left lens, (5.16) 



on the lower lip of the right lens. (5-17) 

Thus S solves the following RH problem: 

• S is analytic on C \ £5, 

• for z G Us, we have S+(z) = S-(z)jg(z), where js is given by 1)5. lljl - Ij5.17|) . 

• as z — > 00, we have S(z) = 1 + Q) . 

Now the ordering of the real parts of the A& in various regions in the complex plane (see 
Lemma l4~2"]) shows that the jump matrices in (|5.13J) - (|5.17|) are all close to the identity matrix 
if n is large, except in a neighborhood of the origin. To be precise, if a < 1 then Re A3 > Re A2 
near the origin in the right half-plane, which means that the entries ± e n ( x 3~ x, 2) i n the jump 
matrices in (J5.14)) and (J5.17|) are not small near the origin but instead grow exponentially if 
n gets large. Similarly the entries ±_ e n ( x ^- x, i) i n the jump matrices in ()5.15|) and ()5.16|) also 
grow exponentially near the origin. On the other hand, if a > 1, then ReAi is bigger than 
the other two in the exceptional neighborhood of the origin, so that the other non-zero off- 
diagonal entries in the jump matrices in (|5.14j) - (|5.17|) grow exponentially in a neighborhood 
of the origin. For a = 1 there are no such exceptions and all jump matrices in (|5.13|) - (j5.17j) 
are close to the identity matrix if n is large. 

When we wrote that certain entries grow exponentially as n gets large, it was understood 
that the value of a / 1 remained fixed. However, eventually we are going to take a = 
1 + 0(n~ 1//2 ) as n — > 00. Then it will turn out that the possible growth of certain entries 
in the jump matrices is confined to a small enough region near the origin, which shrinks 
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sufficiently fast as n — > oo, so that we can still ignore the jumps (|5.1Hj) - (|5.17jl in the next 
step. 

6. Model RH Problem 

We consider the following auxiliary model RH problem: find M : C \ [— z*,z*] — > C 3x3 
such that 

• M is analytic on C \ [—z*, z*}, 

• for x G (—z*,z*) we have M + (x) = M_(x)j'm(^), where 




j M (x)= I -1 0| , for xG (0,0, (6.1) 

and 

j M (x) = I 1 I , for x G {—z*, 0), (6.2) 




• as z — > oo, 

M{z) = I + o(^J . (6.3) 

This RH problem has a solution, see [7J Section 6], that can be explicitly given in terms of 
the mapping functions Wk, k = 1,2, 3, from ()3.1|) and ()3.2j) . The solution takes the form 

(M l {w l {z)) M l {w 2 {z)) Mxiw^zfiX 
M{z) = M 2 ( Wl (z)) M 2 (w 2 (z)) M 2 (w 3 (z)) (6.4) 
\M 3 ( Wl (z)) M 3 (w 2 (z)) M 3 (w 3 (z))J 

where M 1; M 2 , M 3 are the three scalar valued functions 

, r / s w 2 — c 2 , , . . —i w + c , , , , —i w — c . „. 

Mi(tu) = — - M 2 (w) = —= — M 3 (w) = -= — (6.5) 

w\Jw 2 — 3c 2 v2wvw 2 — 3c 2 v2wvw 2 — 3c 2 

Note that by (|3.2|) we have that iuj^z) is of order 2; 1 / 3 as 2 — > 0. By (|6.4|) and (|6.5|) this 
implies that 

M(z) = 0(z- 1/3 ) as z -> 0. (6.6) 

The RH problem for M easily gives that detM(z) = 1. Thus M _1 (2;) exists for z G C \ 
[—2*, z*] and from (j6.6|) it follows that M~ l (z) = 0(z~ 2 / 3 ). However, the special form of the 
solution ()6.4j) - (J6.5|) shows that all cofactors of M are actually 0(z~ 1 ^ 3 ) as z — > 0. Thus 

M-\z) =0(z' 1/3 ) as2->0. (6.7) 

This may also be understood from the fact that M _1 = M*, since together with M it is easy 
to see that also M~* is a solution of the RH problem (|6.1|) - (|6.3j) . 

The model solution M will be used to construct a parametrix for S outside of small 
neighborhoods of the edge points and the origin. Namely, we consider disks of fixed radius 
r around the edge points and a shrinking disk D(0, n~ 1 ^) of radius n~ l / A around the origin. 
At the edge points and at the origin M is not analytic (it is not even bounded) and in the 
disks around the edge points and the origin the parametrix is constructed differently. 
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7. PARAMETRIX AT EDGE POINTS 

The construction of a parametrix P at the edge points ±z* can be done with Airy functions 
in a by now standard way, see jHl E3 HQ E|- We omit details. We only note that ±z* = 
±^^c depends on c and therefore on a. As a — >• 1 we have c — > 1 and so ±z* — > ±^^- We 
construct the Airy parametrices in fixed neighborhoods -D(±^p,r) of so that 

• P is analytic on D(±^fl, r) \ E s , 

• for z G D(±^, r) n S 5 , we have 

P + (z) = P_(s)i5(*), 

where j'5 is given by (|5.11|) - (j5.17|) . 

• as n — > 00, 

P(z) = M(z) (J + O uniformly for 



z± 



3V3 



r. 



(7.2) 



8. Parametrix at the origin 

The main issue is the construction of a parametrix at the origin and this is where the 
Pearcey integrals come in. For a sufficiently close to 1, we want to define Q in a neighborhood 
-D(0, r) of the origin such that 

• Q is analytic on D(0, r) \ 

• for z G D(0, r) fl £5, we have 



Q + (z)=Q„(z) Js (z), (8.1) 

where js is given by ()5.11|) - ([5.17j) . 
• as n — > 00, and with a = 1 + 0(n -1 / 2 ), we have 

Q(z) = M(z) (I + O (n- 1/2 )) uniformly for \z\ = n~ l/i . (8.2) 

The parametrix Q will be constructed with the aid of Pearcey integrals. 

To motivate the construction, we note that the jump matrices for S can be factored as 



3 s 



e- nA -j° s e nA +, 



13) 
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where A = diag(Ai, A2, A3) and 



on(0,**)> (8.4) 



on(-z*,0), (8.5) 



on the upper lip of the right lens, (8.6) 



on the upper lip of the left lens, (8-7) 



on the lower lip of the left lens, (8.8) 



on the lower lip of the right lens. (8-9) 



We show in the next subsection that the Pearcey integrals satisfy a RH problem with 
exactly the above jump matrices except that these jumps are situated on six rays emanating 
from the origin. 



8.1. The Pearcey parametrix. Let b G 1R be fixed. The Pearcey differential equation 
P"'(C) = Cp(C) + bp'iO admits solutions of the form 

Pj(0= [ e-i s4 -^ 2+is<: ds (8.10) 

for j = 0, 1, 2, 3, 4, 5, where 

To = (—00, 00), Ti = (ioo, 0] U [0, 00), 

r 2 = (ioo,0] U [0,-oo), r 3 = (-ioo, 0] U [0,-oo), (8.11) 

T 4 = (—ioo, 0] U [0, 00), T 5 = (—ioo, ioo) 

or any other contours that are homotopic to them as for example given in Figure |H1 The 
formulas (|8.11|) also determine the orientation of the contours Tj. 
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Figure 8. The contours Tj, j = 0,1, ... ,5, equivalent to those in ()8.11|) . that 
are used in the definition of the Pearcey integrals pj. 



Define $ = b) in six sectors by 



~P2 Pi P5 
-P2 Pi PS 
< Pi Pi 

'Po Pi Pa 
Po P'l Pi 
K p'o' Pi Pi 



-P3 



-Pb Pi 



-Ps -Pb Pi 

-p'l -p'l pI 



P4 ~P5 P3 

Pa -Pb Ps 

k p'L -p'l p'l 

% P2 P3 S 

PO P2 PS 

yPg P2 P'l 

'Pi P2 Pb 

P'l P2 Pi 

,P'[ P'i P'l 



for < argC < 7r/4 



for 7r/4 < argC < 37r/4 



for 37r/4 < argC < n 



for — 7r < argC < — 37r/4 



for - 37r/4 < argC < — 7r/4 



for — 7r/4 < argC < 



(8.12) 
(8.13) 
(8.14) 
(8.15) 
(8.16) 
(8.17) 



Then $ has jumps on the six rays. We choose an orientation on these rays so that the rays in 
the right half-plane are oriented from to oo, and the rays in the left half-plane are oriented 
from oo to 0. 
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Then the integral representations (|8.1()j) - (|8.11jl easily imply that $ + = <&-j$ where 

/0 1 0\ 

onargC = 7r/4, j$ = 1—1 on argC = 0, (8.18) 

\° V 
A o o\ 

on argC = -vr/4, >= 1 on arg( = 3vr/4, (8.19) 





on argC = -tt, = 1 on argC = -3vr/4. (8.20) 



So these are indeed the jump matrices of ()8.4|) - (j8.9j) . 

8.2. Asymptotics of Pearcey integrals. A classical steepest descent analysis of the in- 
tegral representations gives the following result for the asymptotic behavior of as 
( — ► oo. As always we use the principal branches of the fractional powers, that is, with a 
branch cut along the negative axis. 

Lemma 8.1. For every fixed b G C, we have as ( — > oo, 



\ -u ^ 1\ /e^(ft») 

$(C;&) = J^ieV 10-11 1 (j + 0(r 2/3 )) e 9 ^^) 

V 3 V C 1/3 / \-^ 2 w 1/ V e^^, 

(8.21) 

/or Im£ > 0, and 

/o- 2 /C 1/3 \ /cu 2 w l\ /e 02 ^'") 

$(C;6) = J— ze^ 10 11 1 (/ + 0(C 2/3 )) 

V 3 V C 1/3 / \^ ^ V \ e e ^ 6 ), 

(8.22)' 

/or Im£ < 0, where to = e 2m ^ 3 and 

e k (C b) = ^ 2fc C 4/3 + ^ fc C 2/3 , k = 1,2, 3. (8.23) 

T/ie O -terms in (|8.21j) and fl8.22j) are uniform for b in a bounded subset of the complex 
plane. 

Proof We give an outline of the proof; cf. also the calculations in jTHj. Let 9(s;(,b) = 
— \s A — |s 2 + i(s. The saddle point equation for ()8.10|) is 

06 

— = -s 3 -bs + i( = 0. 
as 

For b = there are three solutions si = —iu k ( 1 ^ 3 , k = 1,2, 3, and as Q — ► oo, while 6 remains 
bounded, the three saddles Sk = Sfc(C; b) are close to s° k , and in fact 

s k ((; b) = -iu k ( 1/3 - ^C 1/3 + 0(C- 5/3 ) as C - oo. 
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The value at the saddles is 

0(s k ((; 6); £ b) = ^C 4/3 + ^ 2fc C 2/3 + j + 0(C 2/3 ) as C - oo. 

Then, if is the steepest descent path through s fe , we obtain from classical steepest 
descent arguments 



1 «4 b „2 



2tt 



e 4- 2 ->+«°ds = ±J ^ e 0(«*<C*K,6)(i + O ( C -2/3)) 

V-0MC,fe);C,6) 

±^^ 2fc C- 1/3 ef- fc « 4/3+ ^ 2fc C 2/3 +^(l + 0(C- 2/3 )). 



The choice of ± sign depends on the orientation of the steepest descent path. 

Now take any of the six sectors that appear in the definition (|8.12jl - (j8.17jl of $ and take 
some Pj that appears in the definition of $ in that sector. The contour R, in the definition 
(j8.10|) of pj can be deformed to the steepest descent contour through one of the saddles, or 
to the union of two or three such steepest descent contours. However, in the latter case, 
it turns out that there is always a unique dominant saddle for pj in that particular sector. 
Thus for some k and some choice of ± sign, we have 



Pj (0 = ± j|^ 2 rV3 e i^/3 + §^/^ (1 + 0(C - 2 / 3)) (8 . 24) 



as ( — > oo in the chosen sector. Similarly, 

p'jiC) = ±^f ieVC"+V*?">+$ (1 + 0(C 2/3 )), (8.25) 

p"j{Q = ±v /|^ fc C 1/3 e^^ 4/3+ ^^ 2/3+ -(l + 0( C -/ 3 )). (8.26) 

A further analysis reveals which value of k and what sign is associated with pj in the par- 
ticular sector. We will not go through this analysis here, but the result is given by ()8.21j) ) 
and (J8.22|) . This completes the proof of the lemma. □ 

Note that in the above lemma we only state the leading term in a full asymptotic expansion, 
which is enough for the purposes of this paper. We also stay away from situations where 
saddles coalesce. For more asymptotic results on Pearcey integrals in various regimes, see 
[H and the references cited therein. 

8.3. Definition of Q. We are going to define the local parametrix Q in the form 

Q{z) = £(z)$(n 3 / 4 C(,2); n 1 / 2 6(z))e nA(z) e" n22/6 , A = diag(A 1; A 2 , A 3 ), (8.27) 

where E is an analytic prefactor, z i— > ((z) is a conformal map from a neighborhood of in 
the z-plane to a neighborhood of in the C-plane, and z i— > b(z) is analytic. 

We choose ((z) and b(z) so that the exponential factors in the asymptotic behavior of 
$(ri 3 / 4 £(z); n l l 2 b(z)) are cancelled when we multiply them by e nA ^ z ' > e~ nz //6 . We use the 
functions f% and from Lemma 14.11 in the following definition. These functions depend on 
a, and to emphasize the a-dependence we write /^(z; a) and g%{z; a). The functions ((z) and 
b(z) also depend on a. 



22 PAVEL M. BLEHER AND ARNO B.J. KUIJLAARS 

Definition: For z in a sufficiently small neighborhood of 0, we define 

C(^)=C(^;a)=^[/3(^;a)] 3/4 (8.28) 

and 

In (|8.28jl and (|8.29jl the branch of the fractional powers is chosen which is real and positive 
for real values of z near 0. 

Lemma 8.2. (a) There is an r > and a 5 > so that for each a G (1 — 5,1 + 5) 
we have that z h-> ((z; a) is a conformal map on the disk D(0,r) and z i— > b(z; a) is 
analytic on D(0, r). 
(b) In addition we have 

b(z;a) =0(a-l) + 0(z 2 ) as a -> 1 and z -> 0. (8.30) 

Proof. Following the constructions of fj and gj for j = 1,2,3 in Lemmas 13 .2\ 13 .3[ and 14.11 
and their proofs, we easily see that 

h{z;a) = / 3 (z;l) +0(a- 1), g 3 (z; a) = g 3 (z; 1) + 0(a - 1), as a -»• 1, (8.31) 

uniformly for z in a neighborhood of 0, and 

/ 3 (z;l) = l + 0(z 2 ), ^ 1) = 0(z 2 ) as^O. (8.32) 

Both parts of the lemma follow from (l8~3TJ) and and the definitions 1H72H1) and (jOa . 

□ 

From now on we assume that \a — 1| < 5, where 5 > is as in part (a) of Lemma 18.21 
so that z i — ► C(- 2 ! a ) is a conformal map. Near we choose the precise form of the lenses 
so that the lips of the lenses are mapped by z > ((z;a) to the rays arg£ = ±7r/4 and 
arg<^ = ±37r/4. Then from the fact that the jump matrices ()8.18|) - (j8.20|) of $ agree with 
those in ()8.4j) - (j8.9|) . it follows that the jump condition ()8.1|) for Q is satisfied. This holds for 
any choice of analytic prefactor E that is used in (J8.27)) to define Q. We are going to define 
E so that the matching condition (J8.2j) is satisfied as well. 

8.4. Matching condition. To obtain the matching condition (J8.2)) we first note that the 
definitions (J8.28|) and (J8.29|) give us (we drop the a-dependence in the notation) 

C{z) 4/3 = z 4/3 f3{z)j 6 ( z ) C ( z )a/8 = ft ,( z ). 

Hence by (|4.4|) and ()8.23|) we have for Imz > with \z\ < r, 

e k (((z); b(z)) + X k (z) - z 2 /Q = 0, k = 1, 2, 3, (8.33) 
while for Imz < with \z\ < r, 



*a(C(*) 


;b(z)) + X 1 (z) 


- z 2 /6 = 


Ai- 


_(0)=7TZ, 




;b(z)) + X 2 (z) 


- z 2 /6 = 


A 2 - 


.(0) = -iri 




;b(z)) + X 3 (z) 


- z 2 /6 = 


A 3 - 


(o) = o. 



Assume a = 1 + 0(n 1 ^ 2 ). Then it follows from ()8.30|) that 

n 1/2 \b{z;a) \ < C for \z\ < 2^ 1/4 (8.35) 
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for every n large enough, with a value C that is independent of n. As a consequence we can 
use the expansions (|8.21|) . ()8.22|) as n — > oo, because of Lemma |8~T1 We find from ()8.21|) and 
(I8.22|) and the relations ()8.33|) and ()8.34|) between 9% and that the exponential factors in 
the asymptotic behavior (as n — > oo) of 

$(n 3 / 4 C(^);n 1 / 2 6(z))e nA e- n22/6 

cancel if we take z so that 0.9n -1 / 4 < \z\ < l.ln -1 / 4 . So we have proved the following. 

Lemma 8.3. Let a = 1 + 0(n -1 / 2 ). Then we have as n — > oo, uniformly for z so that 
Q.9n~ 1 ' 4 < \z\ < l.ln" 1 / 4 , that 



Q{z) = E(z)$(n 3/A ((z);n 1/2 b(z))e nA(z) e- nz2/6 

'n~ x ^ 




?l ie nb(*r/S E ( z] 








1 1 

n l ' A , 



K(az))(I + 0(n-^)) 



(8.36) 



where 




for lm( > 0, 



for Im( < 0. 



(8.37) 



Proof. This follows from the asymptotic behavior ()8.21|) and ()8.22j) . since we have shown in 
the above that the exponential factors in (|8.21|) and ()8.22|) are cancelled when we multiply 
them by e n H^) e -nz 2 /6_ 

As for the O-term, we note that n 3 / 4 ((z) = 0(n 1 ' 2 ) if 



cn 



- 1 / 4 with 0.9 < c < 1.1, 



so 



that the 0(C" 2/3 ) term in (J8TT^ - (J8~22|) leads to the 0(n~ 1/3 ) term in (fOfij) . □ 
In order to achieve the matching ()8.2|1 of Q(z) with M(z) we now define the prefactor E 



by 



E(z) 




2tt 



ie~ nb{z)2/% M{z)K{C{z)Y l 



n 1 / 4 
1 
n -1 / 4 



(8.38) 



Then the matching condition (J8.2)) follows from (J8.36|) and (|8.38|) . 

It only remains to check that E is analytic in a full neighborhood of the origin. This 
follows since M and K satisfy the same jump relations on the real line. Indeed we have from 
the expressions (|8.37j) for K, for real ( with ( > 0, 




—UJ UJ 2 1" 

-1 1 1 

— UJ 2 UJ 




24 



PAVEL M. BLEHER AND ARNO B.J. KUIJLAARS 



while for real ( < we have to take into account that £ 1//3 and £ have different it- 
boundary values, so that for ( < 0, 





UJ 






(I 


1 




1 '( 




uj 2 








UJ 






(■ 


1 








uJ 2 








-UJ UJ 

-1 f 



-UJ 



UJ 



-1 0, 

These are indeed equal to the jumps satisfied by M; see (j6.1|) and ()6.2|) . Since ((z) is a 
conformal map on D(0,r) that is real and positive for z G (0, r), and real and negative for 
-2 G (— r, 0), we find that M(z)K(((z))~ 1 is analytic across both (0, r) and (— r, 0). Thus 
is analytic in D(0, r) \ {0}. The isolated singularity at is removable, since the entries 
in M(z) and K(((z))~ 1 have at most z _1 / 3 -type singularity at the origin, and they cannot 
combine to form a pole. The conclusion is that E is analytic. 

This completes the construction of the local parametrix Q at the origin. 

9. Final transformation 
We now fix b G R and let a = 1 + 77^. Now we define 



R(z) 



' S{z)M{zy\ for z G C \ S s outside the disks D(0, n" 1 / 4 ) and L>(±^ 

S(z)P(z)-\ forzGD(±^,r)\E S) 
[ S'^)Q(z)- 1 , for ^ G D(0, n-V*) \ £ s . 



2 ' r J' 



(9.1) 

Then is analytic inside the disks and also across the real interval between the disks. Thus 
R is analytic outside the contour shown in Figure M 





Figure 9. The contour The matrix- valued function R is analytic on 
C \ The disk around has radius n -1 / 4 and is shrinking as n — > 00. The 
disks are oriented counterclockwise and the remaining parts of £r are oriented 
from left to right. 



Lemma 9.1. We have R + = R-jn where 

Jr( z ) — I + 0(n _1 ) uniformly for 



3VS 



Jr( z ) = I + 0(n 1 ^ 6 ) uniformly for \z\ = n 1//4 , 



(9.2) 
(9.3) 
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and there exists c > so that 



-en 



2/3 



3r{z) = I + O \ — - uniformly for z on the remaining parts o/£r. (9.4) 

1 + \z\ 

Proof. The behavior (|9.2|) of the jump matrix on the circles around the endpoints is a 

result of the construction of the Airy parametrix. It follows as in [Tfi| IT?] . 

The jump matrix for \z\ = n^ 1 ^ is by (j9.1|) and (|8.2|) (we use positive orientation) 

j R = MQ-1 = M(I + 0(n~ l/3 ))M~ l = M + MO[n~ l/z )M~ 1 {z). 

Since M(z) = 0{z' 1 / 3 ) and M~ 1 (z) = 0{z- x ' z ) as z -> by and (jHUJ, we obtain Q. 

The jump matrix j'_r(z) on the remaining part of S/j is J + 0(e~ cn ) if z G X_r stays at a 
fixed distance of and ±%r^- But now the disk around is shrinking as n increases, and so 
we have to be more careful here. We note that the jump matrix is 

j r (z) = M(z)js(z)M-\z) 

and we want to know its behavior as n — > oo for z on the lips of the lenses near and 
\z\ > n" 1 / 4 . 

The jump matrices js in ()5.14|) - ()5.17j) contain off-diagonal entries ±e n ( Afe ~ A A For a = 1 
these entries are decaying on the contours and so we have for some positive constant c\ > 0. 

Re((A J -A fc )(^;l))>c 1 |z| 4 / 3 

for z on the lips of the lenses near 0. Since Xj(z; a) = Xj(z; 1) + z 2 l 3 0{a — 1) as a — > 1, we 
then get that 

Re ((Aj - X k ){z; a)) > c x z~^ 3 - c 2 \z\ 2/3 \a - 1|. 
Then if a — 1 = (b/2)n~ 1 / 2 and |z| > rir 1 / 4 we easily get that 

Re (Xj - Afc)(z; a)) > c 3 7^ 1/3 

for some positive constant C3 > 0. Then it follows from (J5.14j) - (J5.17|) that 

j s (z) = I + 0(e-^ 2/3 ). 

This leads to (|9.4|) since M(z) = 0(2~ 1 / 3 ) and M~ x {z) = 0(^ 1/3 ) as z -> 0, see (jEHJ) and 
(jSZTJ . □ 

To summarize, we find that i? solves the following RH problem: 

• R is analytic on C \ 

• for z G £_r, we have R + = R-Jr, where jr satisfies ()9.2|) - (j9.4|L 

• as z — > 00, we have R(z) = I + 0(l/z). 

The RH problem for R is posed on a contour that is varying with n. This is a slight 
complication. However we still can guarantee the following behavior of R as n — > 00. 



Proposition 9.2. is n ^ 00 we /iai>e i/jcrf 

1 + u 

uniformly for z£C \ 



/ n -l/6 

i?(z)=/ + - -) (9.5) 
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Since the proof of Proposition 19 .21 is somewhat technical due to the fact that the contours 
are varying with n, we give it in Appendix [X] 

10. Proof of Theorem 11.11 

Now we are ready for the proof of Theorem 11.11 We fix b £ R and take 

b 

a = 1 + 



2Jn 



10.1. The effect of the transformations Y i— ► T i— > S 1 i— ► i?. We are going to follow the 
effect of the transformations on the correlation kernel K n (x,y; a) for real values of x and y 
close to 0. We start from ()2.3|) which gives K n (x,y;a) in terms of the solution of the RH 
problem for Y. The transformation (|5.1|) then implies that 

\n(xl-yi) / e ~n\ 1+ (x)\ 

K n (x,y;a) = ^- -(0 e nX ^ e nX ^) T+\y)T+(x)\ . (10.1) 

2m(x - y) \ I 

According to the transformation T i— > S given in (|5.7|) - (j5.10|) we now have to distinguish 
between x and y being positive or negative. We will do the calculations explicitly for x > 
and y < 0. The other cases are treated in the same way. 

So we assume that x > and y < 0, and both of them are close to 0. The formulas (|5.7jl 
and (|5.9jl applied to (jlO.lj) then give 

in^-y 2 ) / e -nAi+(x)> 

K n {x,y;a) = — (- e nAl +^ e nA3+(y) ) S'+ 1 (y)5 , + (a:) e - nA2 +^ I . (10.2) 



2m(x — y) 







Next we note that for z close to 0, inside the disk or radius n 1 ^ 4 , we have by IjH.ljl . 

= R(z)Q(z) = R(z)Q(z)e nA{z) e- nz2/6 

where 

= Q(2)e-" A(2) e nz2/6 = E(z)^{n 3 ^C{z; a); n^ 2 b{z; a)); (10.3) 
see ()8.36|) . Thus if < x < n -1 / 4 and — rir 1 / 4 < y < 0, we have 

/ e -«Ai+(aO\ /A 

S+(x) I e" nA2 +W I = #(ar)Q(a;) 111 e^ 76 

and 

(_ e nA 1+ (y) o e nA3+(v)) s+\y) = e ny2/6 (-1 l) Q-\y)R- 1 {y). 
Inserting these two relations into (jl0.2)l we find that 

K n( x >y; a ) = lri(*_ V y ) (- 1 X ) Ql\v)Br\y)R(x)Q + {x) flj . (10.4) 

To obtain the scaling limit ()1.6|) of we need the following lemma. 
Lemma 10.1. Let a n = 1 + {b/2)n~ 1 / 2 . 
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(a) Let x n = xn~ 3//4 where xeRis fixed. Then 

lim n 3 ^((x n ; a n ) = x (10.5) 

n^oo 

and 

lim n 1/2 b(x n ; a n ) = b. (10.6) 

n^oo 

(b) Let also y n = yn~ 3 ^ where y G M is fixed. Then 

lim E- x {y n )R- l {y n )R{x n )E{x n ) = I (10.7) 

Proof, (a) Since ((z; a) = z [f 3 (z; a)] 3 ^ 4 by ()8.28|) and f 3 (z; a) — > 1 if z — > and a — > 1 by 
()8.31j) and ()8.32j) . we get that the limit (|l().5jl immediately follows. 

For ()10.6j) we need to go back to the definitions in Lemmas 13.21 13.31 and 14.11 of fi and 
gj, for j = 1,2,3. From I3~2l and its proof it follows that fi(z;a) = /i(0;a) + 0(z 2 ) and 
gi(z;a) = <7i(0) + 0(z 2 ) as 2 -> 0, and the O-terms are uniform with respect to a in a 
neighborhood of 1. Then by (|3.13|) we have 

g 2 (z;a) = g 2 (0;a) + 0{z 2 ) as z -> (10.8) 

uniformly for a in a neighborhood of 1. Since we have; cf. Lemmas 13.31 and (|4.H 



1 /" z 

-z 2/3 g 3 {z;a) = j S - 1/3 g 2 {s]a)ds 



2 

we get from (jl0.8|) that 

g 3 (z;a) = g 3 (0;a) + <D(z 2 ) as z -> (10.9) 

again uniformly for a in a neighborhood of 1. By (|4.5|) we have ^(0; a) = 3c~ 2 / 3 (c 2 — 1) where 
c= (o+V^+8)/4 = (a-l)/3+£((a-l) 2 ) as a -»• 1. Thus c/ 3 (0;a) = 2(a- 1) +0((a- 1) 2 ) 
as a — > 1, and it follows from ()10.9|) and the definitions of a n and x n that 

n 1/2 g 3 {x n ;a n ) = n 1/2 g 3 (0; a n ) + 0(^ 1 ) = 2(a n -l)+0(^ 1/2 ) = 6+0(^ 1/2 ) as n -> oo. 

Then ()10.6|) follows because of the definition ()8.29|) and the fact that f 3 (x n ;a n ) — > 1 as 
n — > oo. 

(b) Since M (z)if (£(z)) _1 is analytic in a neighborhood of the origin and x n — y n = 

0(n -3 / 4 ), we have 

Kidy^Miy^MMKiCixn))- 1 = I + 0{n~ 3 ^) 
as n — > oo. Hence by (18.38)) 

/n" 1 / 4 \ /n 1 / 4 

^- 1 ( 2/n )E(x n ) = e ^ b (y^) 2 -b(x n ,a n f)/8 I j q (/ + 0(n- 3 / 4 )) 1 

\ n 1 / 4 / \ n- 1 / 4 / 

(10.10) 

Note that both nb(y n ; a n ) 2 and nb(x n ; a n ) 2 tend to b 2 as n — > oo because of (jlO.fcij) . Thus we 
also get from (|8.38|) that 

E(x n ) =0(n 1/4 ), and E~ l (y n ) = 0{n l ^ A ). (10.11) 
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Next, we get from ()9.5|) and Cauchy's theorem that for z = 0(n~ 3//4 ) we have 
d ID/ ~\ _ 1 f R ( s ) /-V„-l/6 



, R(z) = / ^^ds = 0(n~ 1/b ) as n -> oo. 

2m J| s | =n -i/ 4 s - 2; 

Then by the mean-value theorem, 

i2(rc n ) - i2(y n ) = 0{(x n - ^n" 1 ' 6 ) = Ofa" 11 ^) 

so that 

R-\y n )R{x n ) =1 + R-\y n ){R{x n ) - R{y n )) = 1 + Ofr" 11 / 12 ). (10.12) 
Combining (fTTTKjl . (fm~TTHl. (fTTHTD and (fTTH2D we obtain ffTTT71 . □ 

Now we can compute the double scaling limit of K n . Indeed, it follows from pO.Hjl . (|1().4)1 . 
and Lemma flO. II that 

lim 4r k » (4i > 4*; 1 + oV) = 6 ) ( 10 - 13 ) 

n^oo rr/ 4 yrr/ 4 rr/ 4 2yn / 

where 

K cusp (x,y;b) = ^ (-1 l) ^(y; 6)$+(x; 6) | 1 | if ar > and y < 0. 

2m(x -y) I q / 

(10.14) 

Similar calculations show that the limit ()10.13|) exists for all x and y, and 



K cusp (x,y;b) = — -(-1 1 0)$; 1 (y;&)$ + (x;&) [ 1 | if x > and y > 0, 

27u(a; —y) 11 



T 
.0, 



(10.15) 



K c " sp (x,y;6) = ^-^-^(-l 1 0)$; 1 (y;6)$ + (x;6) (0| if z < and y > 0, 



(10.16) 



K cusp (x,y;b) = ^ (-1 1) $7%; 6)$ + (ar; 6) I I if 2 < and y < 0. 

2ni{x — y) 11 



(10.17) 



10.2. Different formula for K cusp . To complete the proof of Theorem ll.il we show that the 
formulas (jl0.14j) - (jl0.17|) for K cusp can be rewritten in the form (J1.7)) given in the theorem. 
This involves the Pearcey integrals p(x) and q(y) of p. 8)1 . 
Define 

[Po Pi Pa\ 

$ = Po P'l P4 • (10-18) 
U Pi P47 



LARGE n LIMIT OF GAUSSIAN RANDOM MATRICES WITH EXTERNAL SOURCE. PART III 29 

Then by (18.13)1 we have that $(C) agrees with $(C) in the sector tt/4 < arg£ < 3rr/4, but 
(110.18)) defines <3> in the full complex (-plane, and in particular on the real axis. 

Using the jump relation $ + = for aig( = 7r/4 and arg£ = 37r/4, see (|8.18|) and 

dHUH) , we find that 

~ ^ ° °\ 

$ + (x;6)=$(a;;6) -1 1-1 if x > 

V° i/ 

and 

/I o 0\ 
= 10 if x < 0. 

V-l -1 1/ 

Inserting this into ()10.14|) - (jl0.17)) we find that all four cases lead to 

T 

K cus P(x,y;b) = — -(0 1 l) ^(y; b)$(x; 6) [ I (10.19) 



2iri(x — y) 



.0, 



which is the same expression for all x, y E K. 

Our next task is compute The inverse of $ is built out of solutions of 

= -zg(z) + bq'(z). (10.20) 

It is easy to see that for any solution q of (jl0.20j) and any solution p of the Pearcey equation 

p"\z) = zp(z) + bp'(z) (10.21) 
we have {jpq" — p'q' + p"q — bpq) = so that 

[p, q] := pq" — p'q + p"g — bpq = const. 

It follows that each row of has the form [q" — bq —q' q) for some particular solution 
of (jl0.20j) . More precisely, since <£> is given by ()10.18j) . we have 

(q'{-bq x -q[ gA 
$-1= \q»-bq 2 -q' 2 q 2 \ (10.22) 
\4-bq 3 -q' z q 3 J 

where 

[p , qt] = 1, [pi, gi] = 0, [p 4 , gi] = 0, 

[po, <h] = o, [pi, g 2 ] = i, [p 4 , g 2 ] = o, (10.23) 

[Po, 93] = 0, [pi, g 3 ] = 0, [p 4 , g 3 ] = 0, 
Then if go = g 2 + g3 we have 

[po,go] = 0, [pi,g ] = l, [p 4 ,go] = l, (10.24) 
and from (fiTTTHjl . ffTTTTTT^ . and (fTIT^l) it follows that 

# cusp (s «• b) = Po ^ q °^ ~ Po( x )<io(y) +p'o( x )<io(y) ~ b Po{x)qo(y) , 1Q 25 s 

2iti(x — y) 



30 PAVEL M. BLEHER AND ARNO B.J. KUIJLAARS 

Recall that (jl().21|) has solutions with integral representations 



p{z) 



ds 



(10.26) 



where T is a contour in the complex plane that starts and ends at infinity at one of the 
angles 0, ±7r/2, or tt. Similarly, there are solutions of ()10.20j) with integral representation 



q(z) = 

yv ; 2m 



(10.27) 



where S is a contour in the complex plane that starts and ends at infinity at one of the 
angles ±7r/4 or ±37r/4. 

Lemma 10.2. Let p and q be given by ()10.26|) and ()10.27|) such that T PI £ = 0. Then 

[p, q] = 0. 

J/rnE = {zq} and r and £ intersect transversally at zq, and if the contours are oriented 
so that T meets X in z on the —-side ofE, then 

[p, q] = 1 

Proof. We write [p, g] = pq" — p'q' + p"q — bpq as a double integral, and for convenience we 
take z = 0. So from fTTHfijl and (fiTOTjl . 



\p, q) 



i 

1 

2vri 





_ sf _ s 2 _ h)e ^ +¥ ^L s ^ dsdt 
t 3 + U — S 3 — bs l t 4 + b t 2_l s 4_b s 2 



s Jr 



s - t 



-e 4 



If T PI X = then we can write this as 

1 



2m 




1 1+4,6+2 9 

e 4 2 — 

s — t ds 



e 4 s 2 s 



+ 



2m 




r </£ 



s -t 



- e 4 S 2 S _ 



dt 



g4 t + 2 C 



2 s dsfit 

dsdt 
dtds 



and we can apply integration by parts to both inner integrals. The integrated terms vanish 
because of the choice of contours and the result is 



\p, q] 



2m 
1 

2ixi 




1+4 , b,2_ 1 „4_ b 2 d 
g4 t " r 2 1 ' 4* 2* 



s Jr 



<9s 



1 




1+4 i b+2_ 1 „4_ b 2 (9 
g 4 t T 2 l 4 * 



r Js 



dt 



s-t 
1 

s-t 



dsdt 

dtds = 0. 



If T n X 7^ then we cannot make the splitting of integrals as above, and we have to 
proceed differently. If T and E intersect at zq as in the statement of the second part of the 
lemma, then we can deform contours so that T and £ intersect in 0, and that for some 5 > 0, 
S contains the real interval [—5, 5} oriented from left to right, and T contains the vertical 
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interval [— i8, iS] oriented from bottom to top. Let e G {0,5) and write S e = S \ (— s, e). 
Then it follows as above that 



lim 

e^o 2ni 



lim 



s £ Jr 

-! 

1 



(-t 2 + st-s 2 - b)e 
1 



"2 s dsdt 



r s — t 



-e 4 +2 



ds 



e 4 2* 



dsdt 



1 _I S 4„6 S 2 (9 
g 4 a 2 s 




g4 r + 2 r 


dtds 







If we now do an integration by parts, integrated terms at ±e appear from the second double 
integral. The other terms vanish and the result is 



s + e 



s — e 



1=469 1J, bj 



e 4 ' 



5 £ 



Now we deform T so that instead of the vertical segment [—iS, i5] it contains the semi-circle 
\s\ = 5, Res > 0. Then we pick up a residue contribution from s = e which is equal to 1. 
The remaining integral vanishes in the limit e — > 0, so that we find [p, q] — 1, as claimed by 
the lemma. □ 



Lemma 110.21 allows us to compute $ 1 explicitly. We claim that for j 

1 

2tH 



1,2,3, 



1*4 , 6 t 2_ 

e 4 +2 



-itz 



dt, 



(10.28) 



where Si is a contour in the left half-plane from e~ 3m ^oo to e 37r ^ 4 oo, S2 is a contour in the 
upper half-plane from e m ^ 4 oo to e 3?r / 4 oo, and S3 is a contour in the lower half-plane from 
e~ 37 ™/ 4 oo to e -7n / 4 oo. Indeed, with these contours Sj, and taking note of the definition and 
orientation of To, V±, and T4 in (|8.11jl . we easily get from Lemma 110.21 that the relations 
f)10.23|) hold. Thus for q = q 2 + 53 we find that go = — iq where q is defined as in (II. 8j) . Since 
Po = 2irp, it is then easy to check that the formula ()10.25|) for the kernel is equivalent to the 
formula (|1-7|) in the statement of the theorem. This completes the proof of Theorem 11.11 



Appendix A. Proof of Proposition 19.21 

Let S# be the contour depicted on Figure with orientation from the left to the right 
and in the positive direction on the circles. As usual, we will assume that the minus side of 
the contour is on the right. 

By a simple arc on S# we will mean a connected, 

relatively open, with respect to T, R , subset S^ C T, R , which does not contain any triple 
point of S#, a point where three curves meet. By L 2 (S#) we will mean, as usual, the space 
of measurable functions with 

\\f\\*=(l I/I 2 h) 2 <<*>■ (A.i) 

We have the following general proposition. 
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Proposition A.l. Suppose that a 3 x 3 matrix-valued function v(z), z e Ti R , belongs to 
L 2 (£#) and it is Lipschitz on some simple arc T,° R C S_r. Suppose also that on Y>° R , v(z 
solves the equation 

2tu J TiR z_ - s 

where z_ means the value of the limit of the integral from the minus side, and j R = I + j R 
Then 

R(z)=I-±[ V -^^ds, zeC\Z R , (A.3 
2m Jy„ z — s 



R 



satisfies on Y,° R the jump condition, 

R+(z) = R_(z)j R (z), zeY? R . (A.4 

Proof. From (fA~2|) . (fOjl . 

= u(z), zeSfi. (A.5 
By the jump property of the Cauchy transform, 

R+(z) - R-(z) = v(z)j° R (z) = R-{z)f R {z), (A.6 

hence R+(z) = R^(z)j R (z). Proposition IA. II is proved. □ 
We will solve equation ()A.2|) by the series, 

v(z) = v (z) + v^z) + v 2 (z) + . . . , (A.7 

where 

v (z) = /; Vj (z) = ~ [ V J^M dSj z G Ei?> j > i. (A.8 

We will inductively estimate Vj(z). We begin with some general definitions and results. 
Introduce the operators 

C£v(z) = — - / -^-ds, zeT, (A.9 
2m J r z± — s 

where T is a contour on the complex plane. We assume that v is Lipschitz and L 2 integrable 
if T is unbounded. We have that 

C+-C r =Id (A. 10 

and 

^+ , n- _ n 1 „ - f v ( s ) 



C+ + Cf = C r = v.p. / — ds, zeT. (A.ll 

m J T z — s 

Suppose that the contour Y is given by the parametric equations, 

r = {x = t, y = cp(t), -oo < t < oo}, (A.12 
where (p is uniformly Lipschitz, so that there exists M > such that 

\<p(x)-<p(y)\<M\x-y\. (A.13 
Then as shown in [13], there exists an absolute constant K such that 

||Cr/|| 2 <^o(l + M) 10 ||/|| 2 , (A.14 
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where 

11/112= QV| 2 H) 2 - (A.15) 
This implies similar estimates for C r . If Tq C T then 

C To = PC r P, Pf = X rJ, (A.16) 

hence 

l|Cr ||2< ||C r || 2 . (A.17) 
Therefore, estimate (|A.14J) holds, with the same constant, for any contour 

r = {x = t, y = <p(t), a < t < b}. (A.18) 

Furthermore, it holds, with the same constant, for any complex linear transformation of 
contour (|A.18Jl . Let us denote by Gm the set of all contours which can be obtained by 
a complex linear transformation from a contour (|A.18|) . where (f satisfies (IA. 13j) and is 
different iable. Observe that any interval of a straight line belongs to Go, and any circular 
arc of angular measure less or equal | belongs to G\- 

Suppose now that r = T\ U . . . U T m is a piecewise contour such that 

(1) Tj belongs to Gm, j = 1, • • • , m ', 

(2) the closed contours, Tj and T k , j ^ k, can intersect only at their end-points; 

(3) if Tj and T k intersect then the angle between them at the intersection point is positive, 

Z(r j ,r fc )>e>0. (A.19) 

(4) if Tj and T k , j ^ k, are two infinite contour then they "well diverge" at infinity, so 
that there exists a constant c > such that 

|7i00-7*(*)| >c(H + |t|), (A.20) 

where jj, 7^ are the parametric equations of the contours Tj, Tk, induced by 
parametrization (jA.12|) . 

Theorem A. 2. If T is a piecewise contour which satisfies conditions (l)-(4), then Cr is 
bounded in L 2 , and HCrlh is estimated from above by a constant which depends only on the 
Lipschitz constants Mj of the contours Tj, j = 1, . . . ,m, and on the constants e and c of 
conditions KA.lfy) . I A. 2(J\) . 

Proof. We have to prove that for some K\ > 0, 

|(L7 r /,s)|<^i||/|| 2 N| 2 . (A.21) 

To that end, it is sufficient to prove that for some K 2 > 0, 

|(Cr(xr,/),Xr^)| < K 2 \\f\\ 2 \\g\\ 2 , l<j,k<m. (A.22) 

For j = k, it follows from estimate (jA.13|) applied to a linear transformation of Tj. For j ^ k 
it follows from (IA~T9|) . (IA.2()jl . and the estimate, 

ir r\im^< mMh . (A . 23) 

7T Jo Jo S + t 

Theorem IA. 21 is proved. □ 
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When applied to the contour E_r, Theorem IA.2I gives that there exists a constant K, 
independent of n, such that 

\\Cv B h<K. (A.24) 

By (|A~T0j) . dUU) this implies that 

\\C% a \\2 < K Q = £±± . (A.25) 

From ()A.8|) we have that 

^ = C^Cfe-i)- J>1. (A.26) 

Since 

|&;-i|| 2 < Mdh'-ilh (A.27) 

and 

\\3r\\c = sap \f R (z)\ < Kin - * , (A-28) 
we obtain the recursive estimate, 

\\vjh < Kn-vWvj-iWz, K = K K 1 . (A.29) 

For V\ we have that 

1Mb = \\Cx R U R )h < K \\j R h < K 3 n-^l (A.30) 

Thus, 

\\vjh < K 3 (Kn'^) j n-^. (A.31) 
This implies the convergence of series (IA.7J) in L 2 , for large n. Let us discuss analytic 
properties of the functions Vj. 
Denote 

16 

S« = U S «' (A.32) 
1=1 

the partition of the contour (see Figure EJ) into 16 simple arcs. Fix any e > 0. Let zq be 
any point on T} R such that the distance from zq to the end-points of is bigger than 

e n = erT~±. (A. 33) 

The function j R (z) can be analytically continued from H l R to the ^-neighborhood of the 
point Zo, 

D(z 0} e n ) = {z : dist(z,z ) < e n }. (A.34) 

This implies that 

v 1 {z) = -±- [ M± ds (A.35) 
2m Jz R z^ - s 

can be also analytically continued from H l R to D(zo, £ n ), because we can deform the contour of 
integration, T, R . Then, inductively, we can analytically continue Vj(z) from T, l R to D(z ,e n ), 
by deforming the contour of integration in (jA.8|) . Observe that on the deformed contour we 
have the L 2 -estimate, (jA.31|) . hence by the Cauchy-Schwarz inequality we obtain that 

\vj(z)\ < Kie'^iKn-ty-Vn-*, z G D(z , e„/2). (A.36) 

This proves the convergence of series (jA.7|) in the neighborhood D(zo,e n /2) to an analytic 
v(z). Thus, v(z) is analytic on outside of the triple points. 
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Observe that the function R(z) defined by formula (jA.3|) . denote it for a moment R(z), 
coincides with R(z) defined by ()9.1|) . Indeed, both R(z) and R(z) solve the same RH problem, 
and if Zq is any triple point of E#, then in some neighborhood of zq, 

\R{z)\ < C\z-Zo\~K \R{z)\<C, (A.37) 



for some C > 0. For R(z) it follows from ()A.3|) by the Cauchy-Schwarz inequality, and for 
R(z) it is obvious from ()9.1|) . If we consider now 

X(z) = R(z)R(z)-\ (A.38) 

then X(z) has no jumps on and in a neighborhood of the triple points it satisfies the 
estimate 

\X(z)\<C\z-za\-*. (A.39) 
Therefore, the triple points are removable singularities and X(z) is analytic on C. Also, 
X(oo) = /, hence X(z) = I everywhere on C, and R(z) = R(z). Now we can estimate R(z). 
From fTOfl . 

R{z)=I + y £R j (z), R J (z) = - 7 ^-[ ^fjjMfj ds . (A.40) 

jf=0 ^ 

Suppose that 

dist(z,E fl ) > 0.1n~3. (A.41) 

Then, 

$00 



Z7T 



(is 



< yttt , (A.42) 
1 + z 



and by (TPTJ) . 



K]T\X{Kn 6y +l n b K\{Kn s)% 24 
1 + |z| ~ 1 + |z 

By summing all these inequalities from j = to oo, we obtain that 



_ , « , 111 l' IV l » /to ill IV I » I « ., , 

Rj(z)\ < \ ; ,J = U - ; , J > 1. (A.43) 



?1 6 

= J + O ( ) , dist(.z,£ fl ) > 0.1?i"t. (A.44) 

In fact, the restriction dist(z, £#) > 0.1n~i is not essential, because we can deform the 
contour E#. This completes the proof of Proposition 19.21 
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